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Abstract 

The Dirac equation has been studied in which the Dirac matrices à, (3 have space 
factors, respectively / and /i, dependent on the particle's space coordinates. The 
/ function deforms Heisenberg àlgebra for the coordinates and momenta operators, 
the function f% being treated as a dependence of the particle mass on its position. 
The properties of these functions in the transition to the Schròdinger equation are 
discussed. The exact solution of the Dirac equation for the particle motion in the 
Coulomnb field with a linear dependence of the / function on the distance r to the 
force centre and the inverse dependence on r for the f\ function has been found. 
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Introduction 



The problems with deformed Heisenberg àlgebra with small additions to the canonical 
commutat ional relations have been under a thorough and versatile scrutiny for a period of 
time [1-10]. Deformed commutational relations were studied for the first time in ^1] where 
this issue was raised in connection with the idea of quantisation of space. The question of 
deformation of the Heisenberg àlgebra can be approached along purely practical lines when 
solving eigenvalues problems. When we have a Hamiltonian in the Schròdinger equation 
with the potential not allowing to find the exact analytical solution of the problem we can 
reduce it to a familiar form (for instance to the Hamiltonian of an harmònic oscillator) 
using generalized coordinates and momenta that fail to satisfy the Heisenberg àlgebra. The 
permutation relations between these operators are the so-called deformed relations. With 
this procedure we transfer the "inconvenient" form of the Hamiltonian into a deformation 
of Heisenberg àlgebra. Sometimes this procedure makes it possible to more effectively find 
the approximate solutions of the Schròdinger equation. Some deforming functions allow to 
treat this kind of transfer of inconveniences from the Hamiltonian onto the permutational 
relations as a problem in which the particle mass is position-dependent. 

We can start from the beginning with a "good" Hamiltonian having deformed Heisenberg 
àlgebra with a certain arbitrary deforming function dependent both on the coordinates and 
momentum. Generally speaking, not always will we be able to make the inverse transition, 
i.e. to "toss" this deformation back to the Hamiltonian. For this matter such problems are 
of interest in themselves similarly to those about the motion of a particle with a position- 
dependent mass. At the same, time we have to deal with the problem of mutual ordering 
of the momentum operators and the inverse mass in the kinetic energy. This problem, 
however, does not appear when we resort to the Dirac equation. 

Thus, we arrive at the possibility to formulate the problem about the motion of the 
relativistic particle with a position-dependent mass in the space with deformed Heisenberg 
àlgebra. To study this problem is the aim of this paper. We also give the exact solution 
of the Dirac equation for the motion of a particle in the Coulomb field when its mass and 
deforming function are specifically dependent on the coordinates. 

1 The initial equations 

Let us start from the Dirac equation for a particle with the potential energy U in conven- 
tional notation: 





where à, (3 are the Dirac matrices, the coordinates and momenta satisfy the permutational 
relations with deformed Heiselberg àlgebra: 

Í[xj,Xk\ = 0, 
[xj,Pk\ = ihójkf, n 9 x 
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with the deforming function / = f(x,y,z) dependent on the particle coordimnates only. 
We assume that the particle mass m substituted for a certain effective mass m* is also 
position-dependent: 

m* = mf 1 , fx = fx(x,y,z). (1.3) 

The embedding in the Dirac equation of the functions / and fi implies involvement of extra 
forces acting on the particle alongside of those represented by the function U. We introduce 
a new momentum: 

p = r 1/2 pr 1/2 , 

P = / 1/2 P/ 1/2 , (1-4) 
in the way that coordinates and new momenta become canonically conjugated 

[xj,x k ] = 0, 

[xj,p k ] = ihó jk , (1.5) 
\pj,Pk] = o. 

Now the Dirac equation (jl.ljl looks as follows: 



/ 1/2 (àp)/ 1/2 c + mc 2 /i/3 + f/ ^ = E^. 
We make the transformation 

* = / 1/2 ^, 

as a result of which equation (jl.6j) for the new function will be 



f(àp)c + mc 2 fj + U 



1.6) 



;i.7) 



We can treat this equation as the usual Dirac equation in which the Dirac matrices a. 
are multiplied by certain position-dependent factors: 



ò! = f'ÒL, 



1.9) 



The matrix components òt! and the matrix p are mutually anticommuting. The squares 
of the components of the matrix òt! equal f 2 , and the square of f3' equals f 2 . 

Before we consider the exact solutions of equation (jl.8|) it is expedient to pass to the 
nonrelativistic limit in the Dirac equation in order to find out the properties of the functions 
/ and /i. 
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2 The nonrelativistic limit. The Schròdinger equation 



In order to receive the Schròdinger equation from equation (jl.Hj) at c — ► oo we introduce 
the new function i[) by the following relation: 



4>. 



/(à&c + m^ftp + E-U 

After substituting (|2.1|) in (jl.8j) we find: 

f(àp)f(àp) m 2 c 4 f 2 -{E-U) 2 ihf{àVU) ihcf 



(2.i; 



2m ^ 2mc 2 2mc 

We measure energy from the rest energy mc 2 , 

E' = E-mc 2 , 
and after simple transformations we obtain: 



/(ap)/(àp)_ [J (E'-U) 2 



+ —L(3(àVf 1 )}^ = 0. 



2m 



2mc 2 



mc 



(fï-l)+^fp(àVf 1 )^ = E> 



ihfjàVU) 

2mc 



(2.2) 



From the latter two terms in the parentheses of equation ()2.2j) follows the condition on 
the behaviour of the function f\ in the nonrelativistic limit. Indeed, for the light velocity c 
to drop out of equation ()2.2|) when c — > 00 it is necessary that fi - 1 ~ 1/c 2 . The function 
fi can lead to one at c — > 00 also faster than 1/c 2 leaving no contribution whatsoever in 
the nonrelativistic limit. If 



A 2 -i 



mer 



00, 



(2.3) 



where U\ = U\{x, y, z) is a certain function of the coordinates; then from equation ()2.2|) we 
find its nonrelativistic limit: 



/(q=p)/(q=p) 

2m 



+ U + Ui 



= E'ip. 



We substitute _ 

and assuming that the function / depends on the length r of the radius- vector r after simple 
transformations using the properties of the matrix òc we obtain the following equation: 



1/2" f l/2\2 



2m 



U + AU + uA(p = E'ip 



(2.4) 



mr ar 



(2.5) 



where S = ha/2 is the operator of the particle spin, à = (à x , à y , à z ) are the Pauli matrices, 
L is the angular momentum. 

Expression (|2.4j) can be treated as the Schròdinger equation for a particle with the 
position dependent mass m = m/f 2 where the momentum operator and the inverse mass 
in the kinetic energy operator are specifically ordered: 

i/^ÍpÍ' (2 - 6) 



m^i* -y/m m 

If the particle has a spin then in the nonrelativistic limit the quantity AU remains. We 
refer to it as spin-orbital deformation interaction. 

If we write equation ()2.4|) using the "old" momentum (jl.4j) we have the Schròdinger 
equation in the space with deformed Heisenberg àlgebra: 

p2 \ 

— + U + AU + U 1 )<p = E'cp. (2.7) 

Hence, if in the nonrelativistic theory we start from the Standard Schròdinger equation for 
the study of the behaviour of the particle with the deformed permutative relations (jl.2j) . the 
contribution from the spin-orbital interaction AU gets lost as well as the term U\ caused 
by the dependence of the particle mass on the coordinates (1.3). 

3 The Dirac radial equation 

We consider the particle motion in the central symmetrical field U and functions /, f\ to 
be dependent on the distance r only. We return to the Dirac equation (jl.8j) and reduce it 
to the radial equation. In order to do it we introduce the radial momentum operator. 

p r = r~ (rp — ih) (3-1) 

and a radial component of the matrix òc, 

r 

à r = (cm), n = -. (3.2) 

r 



Further, following [12], we introduce the operator 



hK = p 



a' 



(<t'l) + h 

à 
à 



(3.3) 



and calculating the product à r K we transform equation (jl.8|) into the following one: 

fà r p r c + l —à r pK + mc 2 fj + U J f = £7*. (3.4) 
r / 
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The operator K is the motion integral with the eigenvalues 

k = ±(j + ^j =±1,±2,... , (3.5) 

j is the quantum number of the total angular momentum. That is why in the representation 
where the operator K is diagonal the Dirac radial equation has the form: 

fà r p r c + —tà r f3k + mc 2 /i/3 + U - E J R = 0, (3.6) 
r j 

and 

= YR, (3.7) 

K is the spherical spinor that is the eigenvalue of the operator K, R is the radial function. 
Now we introduce a new radial function R with the following relation: 

R = (^fà r p r c + l -^-à r f3k + mc 2 hÍ3 + E-U^j R. (3.8) 

Substituting this expression into the previous equation ()3.6|) . we find the equation for R: 

(3.9) 

+ ^ C2 { 2fc2 + ihcfà r ^- - ihmc 3 à r Pf^- - {E - Uf\ R = 0. 
r 2 dr dr J 

Here the separation of the space variables from those describing the infernal degrees of 
freedom is possible if 

Cl - (l) = 

dr \ r j dr 1 



al (D = d A 



(3.10) 

dr \r J dr 

where C\, C2 are constants. 

If ()3.1()j) holds then equation ()3.9|) has the form: 

c\fp r f + ftVÀ/£ ( £) + ^Ç^- + m 2 c 4 A 2 -{E-V)*\R = 0, (3.11) 



dr \ r I r 2 

where the operator 



A 7 me 

A = k/3 + —àrd - i—à r f3C 2 . (3.12) 



The operator À does not depend on the radial coordinate and it can easily be reduced 
to the diagonal form. 

Taking into account the properties of the matrices à r from (|3.2|) and (3 we choose a 
representation in which: 









( -i \ 




-? ■ 


à r = 


K i 



(3.13) 



Then operator (|3.12|) is 



A 



k 

\ hc h 



9± + — c ^ 

hc h 



-k 



and its eigenvalues 



(3.14) 



If one works in the representation where the operator À is diagonal, our radial equation 
(|3.11|) finally gets the following form: 



d ff\ h 2 c 2 f 2 k 2 



+ m i cVi -(E-Uy\R = 0. (3.15) 



Let us remark that as the functions /, fi and U are related by two conditions (|3.1U|) . only 
one of them is independent, for instance, it can be the potential energy U. 



4 The Kepler problem 

Now we consider the Kepler problem, that is the motion of the charged particle in the 
Coulomb field when the potential energy 

e 2 

U = —, (4.1) 

r 

where e 2 is the charge square. From equation ()3.10j) we find the deforming function 

/ = 1 + ur, (4.2) 
where v is a constant and the function 

h = 1 + -, (4.3) 
r 

a is a constant and 

Ci = -e 2 , C 2 = a, 
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Then the eigenvalues of À are 



A = ± i / /»- 2 
After a Standard substitution 



/ mca \ '■ 



hc 



where x — x( r )> the radial equation (|H.15j) comes to be 



where 



h 2 d 2 | h 2 | e 

2m dx 2 2mr 2 



dr 

dx = —. 



,„2 



X = E* X , 



From the latter we have 

xv = ln(l + z/r), < x < oc. 
The vàlues with the asterisk in equation (|4.6|) are as follows: 



r ( r + i) = ^ + (=) 



*2 



E" 



E 

mc 2 

E 2 - 



— A - ( 

h 2 k 2 v h 2 v 



fe 2 



\hc 



+ 



m 2m 
m 2 c A h 2 k 2 v 2 



A — mc a, 



2mc 2 

The effective orbital quantum number 



2m 



Vk 2 



a 2 -l 



a 



a — 



Vk 2 ^ 
/mca\ 2 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.Í 



(4.9) 



a = e 2 /hc is the fine structure constant; here the upper value of l* determines the upper 
sign for A (I4.4JI and the lower value sets the lower sign, respectively. Thus, equation ()4.6j) is 
split into two independent equations for the positive and negative vàlues of the quantity A 
from (J4.4j) . If we write the radial coordinate r from equation (|4.7|) explicitly through x and 
substitute r in equation (|4.b|) then after simple transformations we arrive at the following 
equations: 



d 2 A(A - u/2) 



+ 



2B 



dx 2 sinh 2 (xz//2) tanh(xz//2) 



(4.10) 
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where 



2 l*(l* + l) 



A(A-v/2) = v 



B = — + 1/ — ^— : -■ 



2h 2 

2m 



(4.11) 



E* 



h 2 u 2 l*(l* + 1) e*V 



4m 2 

It is well known that this equation has the exact solution ^3] with the energy levels 

2 £2 



V 

A + -rc r 
2 



(4.12) 



{A + vn r /2f 

n r — 0, 1, 2, . . . is the radial quantum number and bound states exist if 

fi > A 2 , A > 0, fi > 0. (4.13) 

As in our case 

A= v -{r + i\ 

then from ()4.12|) taking into account the notations in (|4.8j) for the energy levels fi we find 
the following equation: 

E 2 -m 2 c 4 h 2 u 2 /l2 2 , h 2 v 2 2 
-(k —or) ; — n 



2m 



fi 



— ramc — 



me* 

n = n r + 1* + 1 



m 
h 2 n 2 



4m 

e 2 fi 2 -2N 

— mc a (k + a ) 

mc 2 2m 



n 2 



(4.14) 



is the principal quantum number. 

It is significant that the quantity A drops out of this equation and a dependence on this 
quantity remains only in the effective orbital quantum number l* . Thus, one solution of 
equation ()4.6|) yields the radial function Xn r .i* for l* = \Jk 2 — a 2 — 1 with the energy fi = 
fin.fc! we have the second solution for the negative sign of the quantity À in (|4.4|) and í)4.9|l . it 
equals the function Xn r ,i*+i with the eigenvalue of energy fi n+1;fc . In the nonrelativistic case, 
the first solution gives l* — l — 0, 1, 2, . . ., and the second one l* — l — 1, 2, . . . where / is the 
usual orbital quantum number. Thus the energy levels for the two solutions coincide with 
the exception of the ground state. Here we have the so-called super-symmetry. The Dirac 
equation (jl.8|) for the Coulomb potential with the deforming functions / and fi satisfying 
conditions ()3.1Ü|) reveals supersymmetry. But this issue calls for a separate study. Solving 
equation 1)4.15)1 for fi = E n ^ we finally find 



fi 



ve 



(n 2 + k + a 



n 2 + d 2 



+ 



e 2 a 



/mcy 

\ h ) n 2 + a 2 
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+ 



+ 



mc 



1 + a 2 /n 2 

2 



\ | a 2 | / ve 2 \ 2 / 
n 2 \2mc 2 J \ 



fc 2 + a 2 \/ k 2 + a 2 ' 
1 + au -\ — 



(4.15) 



2mc 



a 

1 + — 



2(k 2 - a 2 ) - n 



(k 2 + a 



^2\2 



1/2 



The condition for the existence of bound states follows from ([4.13)1 : 



Eo V 2 2 

— -e > k + mc a. 

mc 2 m 



(4.16) 



The initial function \l/ contained in equation 1)1.1)1 is found from ()1.7)) . 1)3.7)1 . 1)3.8)1 and 
(ESD: 



$ = /-V2y ( j ór p rC + l J!£à r Í3k + mc 2 h(3 + E-U 



X 



(4.17) 



where x i s the matrix-column with the elements Xn r ,i* and Xn r ,i*+i- 

Formulae 1)4.15)1 — ()4. 17)1 provide the exact solution of the Kepler problem in the Dirac 
theory with Heisenberg àlgebra that is deformed by function (J4.2)) with the position depen- 
dent particle mass in accordance with (J1.3j) . ()4.3)1 . 



5 Discussion of the results 

If in (j4.15j) we put v = 0, i.e. we remove deformation, the energy levels for the Dirac 
charged particle whose mass is position-dependent are obtained: 

mc 2 í me 2 a í fiM 
L= l + a 2 /n 2 [~¥m T+ \ l+ n 2 ~)> (5 ' lj 

and 

e 2 

a < — - 

mc z 

that follows from 1)4.16)1 . This result was originally discovered in [2] and reproduced in 
[15J by a different technique. 

The nonrelativistic limit, c — > oo, for expression ()4.15)) was found. We assume that the 
function fi satisfies condition ()2.3jl . otherwise said, we believe that the dependence of the 
particle mass on its coordinates makes its own contributions into the nonrelativistic limit. 
It means that taking into account the explicit form of the function f\ ()4.3)1 the parameter 
a ~ 1/c 2 . That is why we take 

e 2 

mc z 
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where a is a dimensionless constant. In this case the nonrelativistic limit for the energy E 
is as follows: 

E> = E-mc> = -^(e*-^kA 2 -^+ U -(e> + ^ 

2h 2 n 2 V 2m ) 8m 2\ 2m 

4 

me , , 
+ (5.2) 

and in accordance with ()4.16j) the energy spectrum is limited. 

m 

It is interesting to compare expressions (|5.2j) and (|5.3|) when a = with the results in 
[Tïïj where the Schròdinger equation for a particle in the Coulomb field with the deforming 
function (|4.2|) . in our notation, was solved: 

E' QT = U - — [/(/ + 1) + 1]V - —n* + U -(e* + ^ [1(1 + 1) + l}) , (5.4) 

QT 2h 2 n 2 V 2m L V 7 l J 8m 2 \ 2m L v ; V y 

with the condition that 

e 2 >|^[(/ + l)(2/ + l) + l]. 

The difference of this expression from formula ()5.2|) is explained by the fact that the 
authors of ^B] disregarded the deformational spin-orbital interaction AU which arises nat- 
urally in our treatment in the nonrelativistic limit from the Dirac equation. These authors 
started from the Schròdinger equation at once. If the said interaction is not taken into 
account, we must deduce the contribution from 

AU = —(SL) + —(SL)-, (5.5) 
m m r 

which follows from ()2.5)1 and ()4.2jl . As the eigenvalue of the operator (SL) = (J 2 — L 2 — S 2 )/2 
equals ^ 2 [j(j + l)-/(/ + l)-3/4]/2 = /ï 2 [(j + i/2) 2 -Z(Z + i)-l]/2 = h 2 [k 2 -l(l + l)-l]/2 this 
contribution can be easily taken into account. Indeed, in order to remove the contribution of 
AU from our result it is necessary to deduce from the energy E' the contribution of the first 
term in ()5.5|) equaling h 2 v 2 [k 2 — 1(1 + 1) — l]/2m; the second term in (|5.5J) should be united 
with the Coulomb potential ()4.1j) by the substitution: e 2 — > e 2 + h 2 u[k 2 — 1(1 + 1) — l]/2m. 
Consequently, from (|5.2j) we arrive at expression (|5.4j) . Besides, we must put j = 1 + 1/2 
in condition (J5.3|) limiting the spectrum in expression (J3.5)) for the quantum number k. In 
other words, we should take a higher value of k 2 = (l + l) 2 . 

Now we give the next after the zeroth approximation (j5.2j) term of the development of 
energy by the de grees 1/c 2 . We represent E^ as a sum of three terms: 

= A^W + A2jB (i) + A 3 E {1) . (5.6) 
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The correction does not depend on the parameter v 



L|/f l + a)- i (l + a/3) 



2h 2 n A 

At a = it transforms into the well-known Sommerfeld formula. The correction 

hu \ 2 fi 2 v 2 (n 2 - k 2 Y 



2m 



n 1 



is brought about only by deformation. The cross term 

A 3 £ (1) = [(1 - a 2 )n\k\ - k 2 - nV] + 



+ 



a 2 ) 



2k 4 



-in 
2 V 



n 



H 2 u 2 a 2 
8m n 4 

4 k 4 



y + k 



2\2 



(5.7) 



(5.Í 



(5.9) 



The obtained results are of general interest. They can also be useful for the study of the 
energy spectrum of nanoheterosystems when the electrons mass is position dependent and 
also whenever it is important to take into account relativistic effects, in particular those of 
spin-orbital interaction. 

In the end, let us mention that the question of the application of deformed commu- 
tational relations in Kepler relativistic problem remains open. As non-deformed Kepler 
problem is Lorenz-invariant the question arises whether this property will be preserved in 
the deformed space. Though it is well-known that the quantum space-time with deformed 
Heisenberg àlgebra can be Lorenz-invariant jTTj, it is obvious that in our case the problem 
is not like that. A similar controversy was found in ^7j that studied the Dirac oscillator 
with deformed commutational relations leading to the existance of the minimal length of 
space. However, this problem calls for a more detailed investigation to be suggested in my 
next paper. 

The author is very grateful to V. M. Tkachuk for fruitful discussions. 
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